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ABSTRACT: Considering the non-conservation of the matter content of the universe, a cosmological model has been built 
up, on the basis of an assumption of an inter-conversion between matter and the dark energy. A function of time has been 
introduced to represents the phenomenon of non-conservation of the matter content. By definition, this function is 
proportional to the matter content of the universe at any point of time. Time dependence of this function has been 
determined from Brans-Dicke field equations. The monotonically decreasing nature of this function indicates a continuous 
conversion of matter into dark energy. Time variations of the proportions of matter and dark energy contents of the 
universe have been derived in terms of this function. To find the time evolution of various cosmological quantities, an 
empirical scale factor has been used in the present model. This scale factor has been selected in a way such that it generates 
a deceleration parameter which shows a change of sign with time, implying a change of phase of the expanding universe 
from deceleration to acceleration. Dependence of the scalar field upon the scale factor has also been chosen empirically. 
Time variations of several quantities of cosmological interest have been shown graphically. 
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1. INTRODUCTION 


Several astrophysical observations have confirmed 
that the universe is expanding with acceleration [1 - 3]. 
Studies on cosmology have shown that nearly seventy 
percent of all constituents of the universe have a large 
negative pressure and it is referred to as dark energy, 
which is generally held responsible for driving the 
expansion process with acceleration. 


In general theory of relativity, the cosmological 
constant is a widely known parameter that represents 
dark energy, in spite of its own limitations [4]. The 
researchers have proposed a large number of models 
regarding dark energy and their characteristics have 
been studied in details [4, 5]. The accelerated expansion 
is a very recent phenomenon and it was preceded by a 
phase of decelerated expansion of the universe. This was 
important for the successful nucleosynthesis and also for 
the structure formation of the universe. Astrophysical 
findings show that, beyond a certain value of the redshift 
parameter (z), the expansion of the universe continued 
with deceleration [6]. Only during the later stages of the 
matter dominated era, the evolution of the dark energy 
component of the universe started having its effect on its 
dynamics. Using supernova data, a recent study has 
shown the deceleration parameter (q) changed sign with 
time from positive to negative, indicating the existence of 
a phase of deceleration before the beginning of the 
present phase of accelerated expansion [7]. 


Several theoretical models regarding dark energy have 
been proposed to explain the observational findings [8, 
9]. All these models have successfully generated a 
negative value of the deceleration parameter, implying 
accelerated expansion. There are models involving a 
scalar field with a positive potential which produces an 
effective negative pressure if the kinetic term is smaller 
than the potential term. This scalar field is known as the 
quintessence scalar field. There are a plenty of models 
involving quintessence potentials and they have been 
extensively used. To have detailed information in this 
regard, one may go through a study by V. Sahni on this 
field [10]. A proper physical explanation is often missing 
about the origins of models of most of the quintessence 
potentials. 


It is natural to assume an interaction between the 
main constituents of the universe, i.e. matter and dark 
energy, in the mater dominated era, to explain the 
observations. It leads to a transfer of energy from one 
sector to another. Several models are based on a transfer 
of energy from dark matter to dark energy, in such a way 
that during the later period of evolution, the dark energy 
exceeds matter, causing an accelerated expansion of the 
universe [11, 12]. 


To avoid the difficulties caused by the arbitrariness in 
all these models, non-minimally coupled scalar field 
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theories have been used to effectively explain the change 
from decelerated expansion to accelerated expansion. A 
scalar field is present in its framework. Brans-Dicke 
theory of gravity is the most natural choice in this regard, 
which is a scalar-tensor generalization of general 
relativity. This is mainly due to its simplicity and the 
predictions of results of general relativity in some limit. 
Present acceleration of the universe has been explained 
in the light of Brans-Dicke theory and its modified 
versions [13, 14]. It has been found that Brans-Dicke 
theory can generate sufficient acceleration in the matter 
dominated era without taking into account any 
quintessence field [15]. One actually looks for a theory 
which can account for the transition from the phase of 
deceleration to acceleration. In most of the models, the 
dark energy and dark matter components are allowed to 
evolve independently, without having any interaction. 
Zimdahl and Pavon have shown that the interaction 
between dark energy and dark matter can be used to 
solve the coincidence problem [16]. The idea of an 
interaction or inter-conversion of energy between the 
dark matter and the Brans-Dicke scalar field can 
originate from this concept. The possibility of an inter- 
conversion between matter and the non-minimally 
coupled scalar field was predicted by Amendola [17]. 


Many cosmological models, based on Brans-Dicke 
theory, have shown that the dimensionless parameter w 
must have a very low negative value, of the order of 
unity, for an accelerated expansion of the universe [18]. 
Considering the dark matter to be interacting with the 
Brans-Dicke scalar field, it was once shown that a 
generalized form of the Brans-Dicke theory can account 
for an accelerated expansion even with a high value of w 
[19]. In these models, one either chooses a quintessence 
scalar field or makes a modification of the Brans-Dicke 
theory to account for the findings in the context of the 
accelerated expansion. In a recent study, Clifton and 
Barrow have shown (also shown by Banerjee and Das), 
that no additional potential is required to generate the 
change of sign of the deceleration parameter from 
positive to negative [19, 20]. They considered an 
interaction between the Brans-Dicke scalar field and the 
dark matter to account for the observational findings in 
this regard. 


The present study is based on a generalized form of 
Brans-Dicke theory where the dimensionless parameter 
w is not treated as a constant. It is a function of the scalar 
field g which is a time dependent quantity. This 
generalized form of the theory was first proposed by 
Bergman and it was modified into a form of greater 
applicability by Nordtvedt [21, 22]. This article is not 
based upon any particular theoretical model involving a 
mechanism of interaction between matter and the scalar 
field that causes transfer of energy from one field to 
another. A simple model has been formulated here, only 
by taking into account the fact that the matter content of 
the universe is not conserved. Since the major 
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constituents of the present universe are matter and dark 
energy, the possibilities of an inter-conversion between 
matter and dark energy have been inherently kept open 
in this study. This interaction or inter-conversion must 
have its effect on the parameters that characterize the 
cosmic expansion. The matter content decays into dark 
energy which is regarded as responsible for generating 
the accelerated expansion of the universe, following a 
phase of deceleration. Here we have inserted a function 
of time, denoted by f(t), into the expression 
representing the conservation of matter, to account for 
the non-conservation of the matter content of the 
universe, leading to a new form of the expression for the 
density of matter (9). This new form of the density of 
matter (9) shows that if we choose f(t) = 1 for all values 
of t, it would mean conservation of the matter content of 
the universe. To determine the time dependence of the 
function f(t), we have used empirical expressions of the 
scale factor (a) and the scalar field (g) in Brans-Dicke 
field equations. The function f(t) is defined as the ratio 
of the matter content of the universe at any time (t) to 
the content of matter at the present time. The time 
variations of the proportions of dark energy and matter 
have been obtained in terms of the function f(t), 
assuming the dark energy to be produced from matter. 
The time dependence of dark energy density, total 
energy density and gravitational constant have been 
explored and shown graphically. The main objective of 
formulating the present model is to determine the effect 
of the variation of matter content of the universe on the 
characteristics of cosmic expansion. The consequence of 
generation of the dark energy content upon the time 
dependence of gravitational constant, Hubble parameter 
and the deceleration parameter has been explored in this 
model. The present study shows that the rate of change 
of these parameters increases with the rate of generation 
of dark energy. 


2. Field Equations and the Theoretical Model 


In the generalized Brans-Dicke theory, the action is 
given by [23], 


R 
S=f +O 9, + Lmby—9 atx (1) 


167G 


where R denotes the Ricci scalar, Lm denotes the matter 
Lagrangian, @ stands for the Brans-Dicke scalar field and 
@ stands for the dimensionless coupling parameter 
which is considered to be a function of @ in generalized 
Brans-Dicke theory, unlike its earlier version where it 
was regarded as a constant. 

The line element, for a spatially flat Friedmann- 
Robertson-Walker space-time, is given by, 


ds? = dt? — a*(t) J(r,6,9) (2) 


with J = [dr? +r? do? + r?sin?6 d?] 
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where a(t) stands for the scale factor of the universe. 
The symbols r, 6 and @ denote the spherical polar 
coordinates. 


For a_ spatially flat Friedmann-Robertson-Walker 
space-time (eqn. 2), variation of action (eqn. 1) with 
respect to the metric tensor components yields the 
following field equations of generalized Brans-Dicke 
theory, for a matter dominated (pressureless dust) 
universe [24], 


eet) 3) 
a Neer eae (4) 


Combining equations (3) and (4) we get, 


a ay? _p_ pap @ 
22+4(<) a (5) 


Using equation (3), w(@) can be expressed as, 


oo) =2[3 (2) -2+322] (2) (6) 


ge 2P\\P 


In the above expressions, p denotes the density of 
matter. The content of matter (dark + baryonic) of the 
universe is assumed to remain conserved in several 
theoretical studies [24]. Conservation of the matter 
content of the universe is expressed as, 


pa? = pod) = Po (7) 


We have taken dy = 1 in the above equation. 

In some studies on Brans-Dicke theory of gravitation, one 
takes into consideration an interaction between matter 
and the scalar field [25]. A possibility of an inter- 
conversion between matter and dark energy is taken into 
account in these studies. Here we propose the following 
relation for the density of matter (p), keeping in mind 
this possibility of inter-conversion. 


pa? = f(t) poag = f(t) po (8) 


We have taken dp = 1 in the above equation. 

No theoretical model has been considered in the 
present study to explain or analyze the mechanism of 
interaction between matter and the scalar field. We have 
looked into a simple fact regarding equation (7). Its right 
hand side cannot be independent of time, if one takes 
into account non-conservation of matter due to an inter- 
conversion between matter and dark energy. We have 
incorporated a function of time f(t) into equation (7) to 
get a new relation represented by equation (8). This 


3 
function f(t) = ee is the ratio of the matter content of 
0“%o 
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the universe at any time t to the matter content at the 
present time (t =t )). Thus, the function f(t) is 
proportional to the total content of matter (dark and 
baryonic) of the universe at the instant of time t. We 
propose to denote this ratio by R, where R, = f(t) = 
M(t)/M(t,), with M(t) being a symbol for the matter 
content at time ¢. Let us define a second ratio 

_1df_ 1 dM(t) 

2 fdt M(t) dt 

change of matter content. If, at any instant of time, R> is 
found to be negative, it indicates a loss of matter due to 
its transformation into dark energy. Let us define a third 

‘ M(t)—M(to) 
ratio R3=f(t)-1= ae 
matter content with respect to its value at the present 
time. 


which represents fractional rate of 


indicating a change of 


It has been assumed in the present model that the 
process of conversion of matter into dark energy started 
in the past at the time of t = Atp with A < 1. Therefore, 
M(Ato) = M(to)R, (Ato) was the total matter content of 
the universe, at t = Aty, when no dark energy existed. 
Thus M(Ato) is the total amount of matter and dark 
energy that remains conserved. Assuming the dark 
energy to be created from matter, the proportion of dark 
energy in the universe at any time t is given by the 
following ratio (R,). 


_ MQto)-M(t) _ fAto)-F() 
Ra nts) FA) 0) 


Thus, the percentage of dark energy present in the 
universe at any time is 100 x R,. Approximately 70% of 
the total matter-energy content of the universe is dark 
energy at the present time [18]. For an appropriate 
choice of values of the parameters A and n (to be defined 
later), we must have R,4(to) = 0.7 approximately. 


Since matter and dark energy are the main constituents 
of the present universe, the proportion of matter (dark + 
baryonic) present in the universe can be expressed as, 


Rs =1-R, 
_ _ M(ty)-M(t) __ fe) 
1 NOt) FAte) a 


Thus, the percentage of matter (dark + baryonic) present 
in the universe at any time t is 100 X R(t). 

Our aim is to determine the time dependence of the 
ratios R,,R2,R3,R,and Rs; in terms of the function f(t) 
whose functional form will be determined. 

The density of dark energy, in terms of all these 
parameters, is given by, 


_ Re _ fAte)-FO 
an er C a 


Hence, the density of all constituents (matter and 
energy) is, 
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f(ato) 
Pr=patp=z= a P (11B) 


Using equation (8), f(t) can be expressed as, 
fe) =a & (12) 
Po 


In equation (12), the density of matter (9) can be 
obtained from equation (5). For this purpose, we have 
used an empirical forms of Brans-Dicke scalar field 9, 
following some recent studies in this regard [24, 26]. The 
proposed ansatz for ¢ is given by, 


P= Po (2) = oa” (13) 


In equation (13), n is a constant parameter which 
determines the rapidity with which the scalar field @ 


(= =) changes with time. Larger |n| value means faster 
change in @ with time. We have taken dp = 1 here. 


Using equation (13) in (5), the density of matter of the 
universe (p) is obtained as, 


p = pH? [n? + (4-—q)n+ (4-2q)] (14) 


Using equation (14), the present value of p (denoted by 
Po) can be expressed as, 


Po = PoHo’ [n? + (4 — qo)n + (4 — 240)] (15) 


Substituting equations (14) and (15) into equation (12), 
one gets, 


pH? [n?+(4—-q)n+(4-2q)] 
oHo?[n2+(4-qo)n+(4—-240)1 


f(t) = a* = (16) 

We have used the standard expressions of Hubble 
parameter (H) and deceleration parameter (q), which 
are H=da/a and q=-—da/da? respectively, in our 
derivation of the equation (14). The parameters @, H and 
q in the expression of f(t), in the equation (16), are all 
functions of time. The time evolution of these parameters 
depends upon the time variation of the scale factor from 
which they have to be calculated. We have used an 
empirical scale factor to calculate f(t) using equation 
(16). The criterion used to choose this scale factor is 
based on recent observations regarding the deceleration 
parameter q(=-—da/d?). Universe had a state of 
decelerated expansion before the beginning of the 
present phase of acceleration, as per astrophysical 
observations [18, 19, 24]. Thus, the deceleration 
parameter must have changed its sign from positive to 
negative, in the history of the universe. Our scale factor 
has been so chosen that it shows a signature flip, as a 
function of time. This scale factor, obtained by a little 
variation of the empirical scale factors used previously 
by Roy et al. and Pradhan et al, is given by [26, 27], 
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a = agExpla{(t/to)* — 1}| (17) 


Here a, B are constants which must have the same sign to 
make sure that the scale factor increases with time. In 
terms of this scale factor (eqn. 17), we express the scalar 
field (pg), the Hubble parameter (H) and the deceleration 
parameter (q) in the following way. 


2 = Po (2) 


ao 
= PoExp|na{ (t/t)? — 1}] (18) 
a a, B- 
naten(y 0 
a=- Bett) 20) 


For 0<f6 <1and a>O0, it is found that, gq > +0 as 
t—>Oand q > -—1as t > ©, showing a signature flip of 
q, as per equation (20). Thus, the deceleration 
parameter, obtained from the scale factor (eqn. 17), 
changes sign from positive to negative as time goes on. 
The constant parameters (a,f) have been determined 
from the conditions given below. 


Condition1: H = Hp att = to (21A) 
Condition 2: q=qo att =to (21B) 


Using equations (19) and (20) for the above conditions, 
one obtains, 


Hot 
c= oto 
1—-Hoto(1+40) 


= 1.915 (22A) 


B =1—Hoto(1 + qo) = 0.537 (22B) 


In the present article we have used the following values 
of the cosmological parameters. 


Km 
Hy) =72 (=) per Mega Parsec = 2.33 X 107'8sec1 


ty = 14 billion years = 4.415 x 101’sec 


1 
(3 1.498 x 10°°m-3 Kg s? 
0 
Po = 2.831 x 10-27Kg/m? [present density of matter 


(dark+baryonic)]| 
do = —0.55 


To calculate the value of the function f(t) from 
equation (16), one must use equations (18), (19), (20), 
(22A, B) and also use the above mentioned values of 
cosmological parameters. 


The defining equation of the function f(t) is 
pa? = f(t) poaz. In accordance with this relation, the 
value of f(t)is always positive and, f(t) =1 at t=to 
(taking aj =1). Here, in equation (16), we have 
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f(t) =a)? =1 at t=ty. The ranges of n values for 
which f(t) is positive over the entire range of study (say, 
from t = 0.5tg to t = 1.5t9) are given below. 


n<n, or n > nz where, 
nm, =(q—2)min during 0.5tp to 1.5t (23) 
N2 =(q—-2)max during 0.5to to 1.5t (24) 


Based on the time dependence of q (eqn. 20), during the 
period from 0.5t9 to 1.5to, we find, 


n, = —2.73 and nz = —2.14 


Thus we find that there is a lower and an upper range of 
permissible values for n which aren<n, or n>n, 
respectively. The upper range (i.e. n >n,), has both 
positive and negative values of n and the lower range 
(i.e. n <n,) includes only negative values. According to 
equation (13), the scalar field (@) decreases with time for 
negative values of n, causing the gravitational constant 


(G= 2 to increase with time. Thus it is found that the 


upper range of n allows G to be both increasing and 
decreasing function of time, but G only increases with 
time for the lower range of n values. To choose between 
these two ranges of n, we need to calculate the values of 
Wo for these ranges of n values and compare them with 
those reported by other studies. Using equation (13) in 
(6) we get the following expression of w. 


w =4[3a+n)-S] (25) 


The value of w at the present epoch (i.e. wo), obtained 
from equation (25), is given by, 


2 
Wo ==[34 +n) - 


ae (28) 


Several studies on Brans-Dicke theory show that, wo has 
a negative value close to (—1) [24, 26]. As per equation 
(26), the condition to be satisfied by n, to have wy < 0, is 
given by, 


Po _ 
n< 3boHo? 1 or, n< 0.9884 (27) 


Above condition is satisfied, for the entire lower range 
of n values (ie. n <n,) and for a very small part of its 
upper range (i.e.n > nz). 

Using equation (18), the time varying gravitational 
constant (G), which is the reciprocal of the Brans-Dicke 
scalar field parameter (@), is expressed as, 


G= 


7 5 ExP[—na{(t/to)? — 1}] 


Slr 
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Using the relation Gp = 1/@ we therefore have, 


bei Exp|—na{(t/to)? — 1}] (28) 


a 


A quantity, =, which is obtainable from astrophysical 
observations, is given by, 


6216 nf nH (29) 


According to equation (29) one gets, 


(5) a, = ~"He (30) 


According to Steven Weinberg [28], we must have, 
J] < 4x 10-19 yr, 
Glt=ty 


Hence, the value of the parameter n should be chosen 
accordingly, using equation (30), taking Hp = 2.33 x 
10718sec7}. 


3. Graphical Analysis of Theoretical Findings 


Figure 1 shows the time dependence of the scale factor 
(a) and the redshift parameter (z). Here z = “o = 1; 


where we have used ay = 1 for the entire formulation. 
Several quantities of cosmological interest are often 
plotted as functions of the redshift parameter. Figure 1 
allows one find the value of z at a certain time and also 
the corresponding value of the scale factor. In the figures 
where we have shown the time variations of different 
parameters, we can also find their dependence on z, with 
the help of Figure 1. 


Figure 2 shows the time dependence of the BD 
parameter (w) for three values of the parameter n in its 
lower range. For any value of this parameter, it gradually 
becomes less negative with time. At any fixed instant of 
time, it becomes less negative as we decrease n. This 
behaviour of w is consistent with other recent studies on 
Brans-Dicke theory [29, 30]. 


Figure 3 shows the time variation of the ratio R,(= f) 
for three different values of n in its lower range. For any 
fixed value of n, R, decreases with time at a gradually 
decreasing rate. It decreases at a faster rate for more 
negative values of n. It indicates that a faster rate of 
change of the scalar field (#) causes a more rapid decay 
of matter into dark energy. 


Any increase in the rate of change of @ must be 
associated with a faster change of the gravitational 
constant (G = 1/@). According to equation (30), the 
gravitational constant increases with time for negative 
values of n. Thus, for the entire lower range of n values, 
G increases with time. Only for the negative values of the 
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upper range, which constitute a very small part of this 
range, G increases with time. There are experimental 
observations and theoretical models where G has been 
shown to be increasing with time [27]. The requirement 
of Wo <0 is also fulfilled by negative values of n, as 
shown by equation (27). 

Therefore, it would be logically correct to choose values 
from the lower range of n to determine the time 
dependence of various quantities of cosmological 
interest in the present study. 


Figure 4 shows the time variations of R2 and R3. Here 
Rz is negative over the entire span, implying a decay of 
the matter content into dark energy. It becomes less 
negative with time, implying a decrease in the rate of 
generation of dark energy from matter. R3 decreases 
with time, becoming negative beyond the present time 
(t = to), due to monotonically decreasing nature of the 
matter content with time. 


Figure 5 shows the time variations of the ratios R, and 
Rz. They are found to increase and decrease respectively 
with time. Their rates of change with time become 
smaller as time goes on. Since dark energy is obtained at 
the cost of the matter content of the universe, the sum of 
their proportions, R, and Rs, is unity since the beginning 
of the process of conversion (from matter to dark 
energy) that started at t = Atp. Withn = —4 and A = 0.6, 
these curves give nearly correct values at the present 
time t = to, consistent with recent observations [18]. 


Figure 6 shows the time dependence of the density of 
matter ((,,), the density of dark energy (pg) and the 
density of the total matter-energy content (p,). For 
n = —4and A = 0.6, we get approximately correct values 
of these densities at the present time (t = ty). 


Figure 7 shows the variations of H and q as functions 
of R3. As R3 approaches its present value of zero (from a 
positive value), the dark energy increases with time. The 
rate of its generation from matter decreases with time, as 
shown by Figure 5. It is evident from these curves that 
larger proportion of dark energy causes a decrease in the 
Hubble parameter and the deceleration parameter. 


Figure 8 shows the variations of G/G) and G/G 
functions of R3. As R3 approaches its present value of 
zero (from a positive value), the dark energy increases 
with time. The rate of its generation from matter 
decreases with time, as shown by Figure 5. Figure 8 
shows that as dark energy increases with time, the 
gravitational constant increases. It is found that, as 
|R2| = |f/f| becomes less negative with time, the values 


of G/G decreases. 


Figure 9 shows the time dependence of R,, the 
proportion of dark energy of the universe, for three 
values of the parameter n in its lower range. It shows 
that, for more negative values of n, dark energy increases 
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at a faster rate. It indicates that the higher negative 
values of this parameter cause faster generation of dark 
energy from matter. 


Figure 10 shows the time dependence of Rs, the 
proportion of the matter content of the universe, for 
three values of the parameter n in its lower range. It 
shows that, for more negative values of n, the matter 
content of the universe decreases at a faster rate. It 
indicates that the higher negative values of this 
parameter cause faster decay of the matter content into 
dark energy. 


As n is made more negative, the matter content 
(proportional to R,) decreases at a faster rate (Fig. 3). 
This change causes a faster rise in G (as per eqn. 29). 
Changes of Hubble parameter and the deceleration 
parameter are also connected to the change of R, with 
time. It can therefore be concluded from these figures 
that the dark energy content and its rate of production 
from matter are found to play an important role in 
determining the time dependence of various 
cosmological parameters connected to the expansion of 
universe. 


We carried out a similar research regarding the decay 
of the matter content of the universe into dark energy 
[31]. In that paper we had graphically shown the 
dependence of wy upon the parameter n. That graph had 
a discontinuity at n=0 which is consistent with 
equation (26) of the present article. According to 
equation (13), the value of n = 0 is inconsistent with the 
fact that the scalar field @ is a function of time in Brans- 
Dicke theory. Hence, the value of wy at n=0 has no 
theoretical validity. That article had the plots of Rz and 
R3 as functions of R,, from which their time dependence 
was not very clear. Figure 4 of the present article shows 
their plots as functions of time. In that study we had 
shown the plots of H,q,G/Gp and G/G as functions of R,. 
The present study show their variations as functions of 
R3, in Figures 7 and 8. Here, R3(t) is a function of time 
such that 100 x R3(t) is the percentage of change of the 
matter content at any time ¢, while R,(t) is the 
proportional to matter content of the universe at that 
instant. Time dependence of the Brans-Dicke parameter 
w has been shown in this study in a better way than that 
of our previous study. In Figure 2, it has been shown for 
three different values of the parameter n in its lower 
range, having greater theoretical validity. The present 
study has two additional figures (Figs. 9, 10) showing the 
time evolutions of R, and R; for three different values of 
the parameter n, indicating a faster transformation of 
matter into dark energy for more negative values of this 
parameter. Figure 1 has the plots of time variations of 
the scale factor and the redshift parameter (z). This 
figure enables one to estimate the values of other 
parameters (w, R,,R2,R3,R4, Rs ) for different values of z, 
which was not possible from our previous theoretical 
formulation in this regard [31]. Time dependence of all 
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cosmological quantities discussed here, have been 
determined by using an empirical scale factor (eqn. 17), 
which is different from the one used in our previous 
model. Expressions of @, H and q (eqns. 18, 19, 20),based 
on this scale factor are also different from those obtained 
in the previous study. For both models, the scale factors 
have been so chosen, that the deceleration parameter, 
obtained from it, shows a signature flip from positive to 
negative, implying a change of phase of the expanding 
universe from deceleration to acceleration. Similarity of 
results, between these two models, regarding the time 
variations of the matter content and the dark energy of 
the universe, serves as a confirmation of the basic 
theoretical premise and the underlying assumptions that 
have led to these theoretical formulations, based on the 
field equations of the Brans-Dicke theory of gravitation. 
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Fig. 1: Time dependence of the scale factor (a) 
and the redshift parameter (z), plotted along the 
left and right vertical axes respectively 
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Fig. 2: Time dependence of the Brans-Dicke 
coupling parameter (w) for three values of the 
parameter n 
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Fig. 3: Time dependence of R,[= f(t)] for three 
values of the parameter n 
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Fig. 4: Time dependence of Rz and R3, plotted 
along the left and right vertical axes respectively. 
For this plot, n = —4 


Fig. 5: Time dependence of R, and R; 
(proportions of dark energy and matter 
respectively). For this plot, n = —4 and A = 0.6 
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Fig. 6: Time dependence of the density of matter 


(Pm = p), dark energy (pq) and the total matter-energy 


density (p;). For this plot, n = —4 and A = 0.6 
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Fig. 7: Variations of H and q as functions of R3, shown 


along the left and right vertical axes respectively 
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Fig. 8: Variations of G/Gp and G/G as functions of 
R3, shown along the left and right vertical axes 
respectively. For this plot, n = —4 
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Fig. 9: Time dependence of R, (proportion of 
the dark energy of the universe) for three 
values of the parameter n. For this plot, 2 = 0.6 


Fig. 10: Time dependence of R; (proportion of 
the matter content of the universe) for three 
values of the parameter n. For this plot, 2 = 0.6 


4, CONCLUSIONS 


The aim of the present article is to determine the time 
dependence of the matter content and the dark energy of 
the universe by using standard empirical expressions of 
the scale factor and the scalar field in the field equations 
of Brans-Dicke theory, of a homogeneous and isotropic 
universe with zero spatial curvature. A new function f(t) 
has been assumed to exist in the expression of the 
density of matter. This function is proportional to the 
matter content of the universe. Initially, no particular 
form of this function was assumed. The functional form 
of f(t) has been determined from Brans-Dicke field 
equations. This function has been found to decrease with 
time monotonically, indicating a continuous conversion 
of matter into dark energy. Through a graphical 
depiction of theoretical findings, it has been shown that, 
with the rise of the proportion of dark energy content of 
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the universe with time, the deceleration parameter 
becomes more negative and the gravitational constant 
increases. All these findings of this model clearly indicate 
that the time evolution of the deceleration parameter 
and the gravitational constant depends on the content of 
dark energy and its rate of generation from the matter 
content of the universe. Signature flip of the deceleration 
parameter, from positive to negative, indicates a 
transition of the universe from a phase of decelerated 
expansion to that of an accelerated expansion. This 
change of phase must be connected to the generation of 
dark energy from matter. Considering matter to be the 
only source of dark energy in the present model, it has 
been shown that the present proportions of both of them 
must depend on the duration over which the process of 
conversion (from matter to dark energy) has taken place. 
The smaller the value of the parameter A, longer would 
be the duration of matter-energy conversion. The 
proportions of the matter and dark energy contents at 
the present time (t = to), obtained from the present 
study by a proper tuning of the parameters n and A, are 
consistent with the values obtained from various 
astrophysical observations. As the parameter n is made 
more negative, the function f(t) (i.e. R,) will have a 
faster rate of change with time (as shown in Fig. 3) and 
consequently, R, and R; will have more rapid changes 
with time (in accordance with eqns. 9 and 10 
respectively). It shows that the rate of change of the 
scalar field or gravitational constant is clearly associated 
with the time evolution of the proportions of matter and 
dark energy. One may think of a future plan to improve 
this model by assuming a new empirical expression for 
the Brans-Dicke scalar field parameter (g) and also a 
new empirical form of the scale factor. One may also 
think of defining f(t) in terms of @ and the scale factor. 
The functional form of f(t) should be so derived that its 
value is unity at the present time (t = tg), according the 
expression of equation (8) which serves as a defining 
relation for f(t). Using different functional forms of f(t), 
one needs to find out the nature of dependence of the 
cosmological quantities like Hubble parameter, 
deceleration parameter, gravitational constant etc. upon 
the proportion of matter and dark energy present in the 
universe and also upon the rate of their inter-conversion. 
The present theoretical model provides us with a very 
simple method to determine the nature of time 
variations of the main constituents of the universe 
(matter and dark energy) and also the nature of 
dependence of the time evolutions of various 
cosmological parameters upon the present quantities of 
matter and dark energy in the universe and the rate of 
their inter-conversion. 
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